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I. INTRODUCTION 



The Gutzwiller variational method has been used since long for the study of ground 
state properties of Hubbard-type models for correlated fermions.§ The one-band (one or- 
bital) formulation has been extended in various ways to treat, e.g. antiferromagnetismlH or 
to include additional, uncorrelated orbitals.@~i The Gutzwiller ansatz has also been inves- 
tigated analytically.! However, attempts to extend the method to n = 2 correlated orbitals 
per site have, sofar, led to unsatisfactory results.! 

In this paper we generalize the Gutzwiller method to cases of arbitrary numbers n of corre- 
lated orbitals per site. The extension makes it possible to study the ground state properties 
of realistic models for magnetic 3d transition metals and compounds. As a first applica- 
tion, we examine a two-band model with more general on-site interactions and focus on the 
metal-insulator transition! at quarter and half band fillings. 

II. GENERALIZATION OF THE GUTZWILLER METHOD 

For one band, the Gutzwiller variational wavefunction reads as 
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where j^o) is an uncorrelated wavefunction on a lattice of L sites and the operators D s = 
n s ^n s i measure the double occupancy of sites s. The variational parameter g will reduce 
the total number D = M 12 of doubly occupied sites and may vary between g = 1 (M 12 = 
M\MijV) and g = (M12 = 0). Here Mx^M?) are the gross numbers of spin up(down) 
sites. The norm of (|l|) is approximated by0 

(* I *> = P^Y^g^NniL.M^M^ , (2) 
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with the probabilities Pk = M^.iL — Mf.)\/L\ and the combinatorial weight factors 



Here we have introduced net numbers of occupied sites Mk = Mk — Mi 2 and of empty sites 
Mq = (jj — Mi — M 2 — M12). In the thermodynamic limit the sum of eq. (|2]) is replaced 
by its largest term, leading to the relation 



g 2 = m mi 2 /(mim 2 



with m$ = Mq/L, etc. Also the expectation values 
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are found by combinatorics. Here the "loss factors' 

1 



qk = 



y/m k m + ^m kl mi 
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m k (l - m k ) 

with k, I — 1,2 and k 7^ I represent the step sizes of the momentum distribution functions 
at the Fermi level (qk — 1 for g — 1). 

When we extend the Gutzwiller method to cases of n > 2 correlated orbitals per site, the 
number of different site occupancies J n = 4 n will increase enormously. For n = 2 and using 
the notation af= 1, a|=2, 6|=3, 6|=4 the 16 occupancies are: i) empty (M ); ii) four single 
(M k ,k = 1, 4); hi) six double (M k i ,k,l — l,...,4;k 7^ Z); iv) four triple (M k i p ) and v) 
quadruple (Mi 234 ). In general, the J n occupancies are labeled by subscripts Ij, which are 
assigned to the multiple orbital-spin subscripts by 



{h, I a , J 2n+1 , J 2n+2 , J Jn } = {(0), (1), (2 • n), (12), (k...p)} . 



The net and gross occupancies are measured by appropriate operators, e.g., 
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There are now K n = J n — (2n + 1) multiple occupancies so that the generalized Gutzwiller 
operator reads as 



It fef 



(9) 



j=2n+2 



Note, that in contrast to ref. §, we found it essential to include all multiple occupancies and 
to use net multiple occupancy operators. The norm is given by 



(*i*> = nft£*s n (^p 

k=l g j=2n+2 
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with the P k factors of eq. @, and 



N 5 = L\ 



(11) 



The sum of eq. ( |10D includes all sets of multiple occupancy configurations D = {MjA with 
j = (2n + 2) , J n . Replacing the sum by its largest term leads to the K n relations 
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Keeping the gross numbers m k fixed, we get 



g\ = m m k i/ {m k mi) 



#L P =( m o) 1} m k .. P / (m k ...m p ) , (13) 

where 7 is the number of subscripts in m k .. q . For each of the generalized "loss factors" q k i 
the evaluation of 4 2n_1 terms is required. Nevertheless, the resulting equations for the q k i 
factors look remarkably simple: 

qkk = 77^ r [y/m k m + V ' \Jm k irni 

m k (l-m k ) y 1 



+ \J m ki P m P + "' \l m ki Pq mi Pq + •••] 2 (14a) 

l,p l,p,q 

Qki = VkkQu (14b) 

Here the summation primes exclude k — I, etc. Obviously eqs. ( |T4D are a straightforward 
generalization of eq. (|^) and of results obtained in ref. £|. Simpler forms of the q kk equations 
are obtained, when certain multiple occupancies are put zero, as was assumed in ref. |H| . 



Recently, eqs. (|lj) have also been found by T. OkabeJ 
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Using the generalized Gutzwiller wavefunction, we can now investigate extensions of the 
Hubbard model for arbitrary numbers of orbitals a, (3: 

+ D , J^a+&5^ ) (15) 

ce,/3,s,cr,cr' 

which include, apart from general hopping terms, the on-site interaction U, both orbital- 
diagonal (U aa ) and off-diagonal (ll af3 , a ^ pfj , and the on-site or Hund's rule exchange 
ja/3 r-Q^ p r j mes exclude non-physical orbital-spin combinations in the respective U and J 
summations. For the expectation values of the on-site exchange only the terms 

" Ys'J^KA (16) 

survive, i.e., all spin fluctuations are ignored. 

By assuming J af3 < U af3 < U aa , it is possible to achieve the filling of incomplete d shells 
in the atomic limit according to Hund's rules. Thus, the Hamiltonian of eq. ([15]) represents 
realistic models for the valence electron structure of both 3d transition metal elements and 
of magnetic 3d oxides and halides, provided that a sufficient number of orbitals and sites are 
included. The N— particle wavefunctions , |\l/o) can be chosen to represent para-, ferro-, 
or antiferromagnetic states, and in principle, even more complicated cases, such as orbital 
ordering. Also, ligand orbitals or s, p type orbitals on the d sites may be incorporated. 
However, any additional complexity of the wavefunctions |\l/o) leads to further variational 
parameters in \^>q) representing effective orbital energy shifts or effective hoppings, etc.0 
The problem of orbital density depletion caused by the reduction of multiple occupancies 
can be treated in analogy to ref. |j. 

III. TWO-BAND MODEL WITH ORBITAL DEGENERACY 

In the following, we investigate the n = 2 Hamiltonian, i.e., we assume one atomic site on 
a simple cubic lattice with two, degenerate d(e g ) orbitals, and we study only paramagnetic 
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solutions. The single particle bands are constructed using realistic INN and 2NN hopping 
matrix elements T dda (lNN) = leV, T dda (2NN) = 0.25eV, and T dd5 : T dd7T : T dda = 0.1 : 
—0.3 : 1. There are three interaction parameters U aa = U, U al3 = U' and J af3 = J, which, 
in the limit of vanishing configuration interaction are related by0 

J=\{U-U'). (17) 

For the paramagnetic case we have m k = m, rUup = t, m.1234 = /, m"i2 = m 34 = d d , rriu = 
m 2 3 = d Q , and m 13 = m 2 4 = d t . All g fci are equal 

q = —— r[(\fdd~+ J~d + Jd t )(Vm + Vt) 

mil — m) v v v 



+y/mjm+ ^tf} 2 . (18) 

The energy function for e/a = 4m (number of electrons per atom) is 

£ = 2ge(m) + 2[/d d + 2[/'d + 2(U' - J)d t 

+ (2C/ + 4C/'-2J)(2t + /) . (19) 

Here 2e(m) is the kinetic energy of the uncorrelated case. E has been minimized numerically 
with respect to the five variational parameters d d , d a , d t , t, f and has been studied at several m 
values and for various sets of interaction parameters U and u' = U'/U, keeping 2 J — U — U' 
according to eq. ([17|) . 

There are two m regions of interest. The first is around m = 0.25 (or m = 0.75), the 
other around m = 0.5. For these integer fillings (e/a = 1,2,3) we observe metal-insulator 
(MI) transitions in the (U, u')— plane. However, there are big differences between the two 
cases, mainly concerning the shapes of the phase diagrams in the (U, u')— plane and the 
order of the MI transitions. 

For m = 0.25, the MI phase diagram is shown in Fig. |l]a. For u' = 1, we find Uq ~ 13.8 in 
agreement with ref. [D]. For it' > 1/3 (W > J) the system will avoid all multiple occupancies 
in the limit U — > 00, so that there is always a MI transition. This transition is of second 
order in q, very similar to that of the one-band model.^ Near any critical values (Uc,u' c ), 



q behaves like q ps qo(l — U/Uc), when we approach Uc from the metallic side, i.e.; the 
effective mass m* = m e /q diverges with (1 — U /Uc)^ 1 ^^ For v! < 1/3 (U' < J) the system 
is always metallic. 

In contrast, the MI transition for m = 0.5 is very different (see Fig. |l|a). The transition is 
generally of first order in q. This means that the q values change discontinuously at critical 
values (Uc,u' c ) from a metal (typically q ps 0.6) to an insulator (q = 0) (see Fig. |l]b). Only 
at the point u' = 1 (i.e. Uc = U' c ) the MI transition is second order in q. Note that apart 
from the point u' = 1 the ground state energy E varies at the transition in such a way that 
its first derivatives with respect to all relevant system parameters changes discontinuously. 
In particular we observe a finite discontinuity of the chemical potential /i(m) = dE(m) / 'dm 
at m = 0.5 for all vl < 1 if U exceed Uc by any infinitesimal amount. 

In the insulating region we have t = f = 0, and because of particle hole symmetry, 
m = m = 0. For u' < 1 only the spin triplet occupancy d t prevails: d t = 0.5, d<i = d Q = 0. 
The case u' — 1 is special, as there is no preference for one of the three double occupancies. 
We believe, because of this kind of frustration the metallic phase extends so far out near 
v! < I. 

In the metallic regime all multiple occupancies are nonzero (see Fig. 0). For u' < 1, dt 
increases with U, typically up to d t ~ 0.25. Then, at Uc, d t jumps up to 0.5, representing a 
low-to-high-spin transition. For certain band filling values close to m — 0.5, the first order 
step is still present in the q(U,u f ) and d t (U,u') curves. Further away from m = 0.5, the 
curves smoothen, yet still vary rapidly in the range of critical parameters (Uc, u' c ) as shown 
in Fig. 0. A typical phase diagram in the (U, u')— plane of this "good-to-bad" -metal or 
low-to-high-spin transition is shown in Fig. |[ 

Note that for our evaluation a homogenous phase was assumed. However, in the vicinity 
of the low-to-high-spin transition, the homogeneous phase is unstable against phase sepa- 
ration. As a consequence there exists a surface m ps (U,u') in the space of (U,u',m) which 
divides the homogeneous from the separated region (see Fig. |3|). At any point (U,u',m) 
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with m ps < m < 0.5 (or m' > m > 0.5), the system breaks up into one part with density 
mi = 0.5 and another one with m 2 = m ps (m'). It is an open question, however, whether 
the phase separation is an artefact of Hubbard-type models, because of, e.g., the neglect of 
long range Coulomb interaction. 

The first order behavior of the MI transition for m = 0.5 is caused by the competition 
of two different energy minima, now possible due to of the much larger space of variational 
parameters. We expect that the first order nature will survive when we include antifer- 
romagnetic states in the minimization. First order, low-to-high spin transitions are also 
observed for two-band, itinerant ferromagnetic states.0 In fact, we suspect that first order 
changes in the electronic structure of incomplete 3d shells may drive many phase transitions 
of magnetic 3d compounds. It should be interesting to reexamine the MI phase transitions 
observed e.g. for Ni S, Ni J 2 or RE Ni oJ^t^j\ with respect to their first order behavior. 

Very recently, a dynamical mean field study about the MI transition in a degenerate 
(n=2) Hubbard model has been presented.^ The authors only consider the cases J = 
and/or m < 0.25 (in our notation). In aggreement with our results, they find a qualitative 
feature, which is very similar to the n — 1 case. Nevertheless, their more general remark 
that the nature of the Mott transition is not qualitatively changed by orbital degeneracy is 
not supported by our results, concerning the cases J > 0, m ps 0.5. 

IV. CONCLUDING REMARKS 

In conclusion, we have extended the Gutzwiller method to arbitrary numbers n > 2 of 
correlated orbitals per site and have reported studies on the simplest n = 2 model. We find 
a wealth of new features, not present in the n = 1 models, such as first order metal-insulator 
transitions or low-to-high spin transitions. Numerical studies appear to be feasible even for 
n = 5 realistic models of 3d transition metals and compounds. 
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FIGURES 

FIG. 1. a) Phase diagrams in the (U, u')— plane of metallic and insulating phases for m = 0.25 
(dashed line, LH scale U\) and for m = 0.5 (solid line, RH scale L^). The lines separate the metallic 
phases (small U) from the insulating phases (large U). For m = 0.5 the insulating phase is a triplet 
state (dt = 0.5) for u' < 1, and a singlet state {da = 0.5) for u' > 1. b) Loss factor qc at the MI 
transition points (Uc,u' c ). 

FIG. 2. Double (dt,dd,d ), triple (t) and quadruple (/) occupancies for m = 0.5 (e/a = 2) as 
a function of U, with u' = 0.8 (solid lines). Also shown are dt{U) curves for m = 0.49, 0.48, 0.47, 
respectively, indicating the "good-to-bad" metal transition away from m = 0.5. 

FIG. 3. Phase diagrams in the (U, m)— plane, with v! = 0.8: Lines dividing the "good metal" 
from the "bad metal" phases are indicated by nigh, lines dividing the homogeneous from the phase 
separated regions are denoted by m ps . For the part of m g b given by dashed lines, the discontinuity 
in q(U, v!) has disappeared; here inflection points are used to determine the phase boundary. 
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